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The quantum theory of a harmonic oscillator with a time dependent frequency arises in several 
important physical problems, especially in the study of quantum field theory in an external back- 
ground. While the mathematics of this system is straightforward, several conceptual issues arise in 
such a study. We present a general formalism to address some of the conceptual issues like the emer- 
gence of classicality, definition of particle content, back reaction etc. In particular, wo parametrize 
the wave function in terms of a complex number (which we call excitation parameter) and express 
all physically relevant quantities in terms it. Many of the notions - like those of particle number 
density, effective Lagrangian etc., which are usually defined using asymptotic in-out states - are 
generalized as time-dependent concepts and we show that these generalized definitions lead to use- 
ful and reasonable results. Having developed the general formalism we apply it to several examples. 
Exact analytic expressions are found for a particular toy model and approximate analytic solutions 
arc obtained in the extreme cases of adiabatic and highly non-adiabatic evolution. We then work out 
the exact results numerically for a variety of models and compare them with the analytic results and 
approximations. The formalism is useful in addressing the question of emergence of classicality of 
the quantum state, its relation to particle production and to clarify several conceptual issues related 
to this. In Paper II (arXiv:0708.1237) which is a sequel to this, the formalism will be applied to 
analyze the corresponding issues in the context of quantum field theory in background cosmological 
models and electric fields. 

PACS numbers: 



I. INTRODUCTION AND MOTIVATION 

The quantum harmonic oscillator with a timc-dcp(^nd(mt frequency has been widely studied in numerous contexts, 
because it rears its head in practically every situation involving the study of a quantum field in a non-trivial time 
dependent background (for a review, see for e.g. [1, 2]). The dynamics of such a system is encoded in a propagator 
Kernel K{q2,t2; qi,ti) which will allow us to determine the quantum state ^(92,^2) a,t t = t2 if the quantum state 
tp{qi,ti) a,t t = ti is given. Since the Lagrangian of the system is quadratic in q, the Kernel can be expressed in 
terms of the solutions to the classical equation of motion (see for e.g. [3]). So, if the classical solution for a harmonic 
oscillator with a particular time dependence is known, the quantum theory should be trivial. 

And so one would have thought. The reason one continues to investigate this problem is not so much because 
the equations cannot be solved but because the interpretation of the solutions is non-trivial. (In this paper, for 
example, the solution appears within the first ten equations and the rest of the paper is on interpreting it!) Hence, 
it is worthwhile to raise some questions which need to be addressed in reasonably precise terms, before we proceed 
further. This will provide the motivation for contributing yet another paper on this topic to the already extensive 
literature! 

(1) In virtually every context we are interested in, there will be a classical degree of freedom symbolically denoted 

by a variable C (cosmological background, external electromagnetic field, ) interacting with a quantized degree of 

freedom q (usually a scalar field) with the total Lagrangian for the combined system written as L{q,C) — Lq{C) + 
Lj (C, q) . We will be interested in quantizing q in the background provided by C and study the effect of g on C at the 
semiclassical level. To the lowest order, the configuration of C will be determined by the equations of motion arising 
from Lq, ignoring q completely. If this configuration is nontrivial (say, time dependent) then the quantum theory of 
q will be based on a time dependent Hamiltonian and the q- particles will be generated by the interaction. The first 
conceptual question is: how can one define a notion of such particles and their production when the time dependence 
is nontrivial? If the classical system reaches a time-independent state asymptotically, it is straightforward to define 
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the notion of particles in the asymptotic in and out states and also obtain an expression for the total number of 
particles produced. But many physically relevant problems (like cosmological particle production) will not give us 
this luxury of well defined in and out states. 

(2) The question of defining the notion of a particle in a time dependent background is not one of idle curiosity. If 
C is producing the quanta of q, it has to supply the energy for the process and obviously this will modify the evolution 
of C. It will be important to obtain the equations of motion for C with this backreaction incorporated as the particle 
production progresses. This is certainly not an asymptotic notion and one would like to imagine that - in any causal 
theory - the backreaction on C at time t = ti should not depend on how the system will evolve at time t>t\. So, 
purely conceptually, we need a notion of back reaction which does not use the concept of an out vacuum state etc. 

(3) Such a back reaction — and in fact, the notion of semiclassical evolution — can be meaningful only if we have 
at least an approximate notion of 'classicality' for the q mode. Our intuitive idea of a particle which is produced (that 
has drained away energy from C) is classical and one assumes that a particle which is produced 'stays produced'. 
This is of course not true in general and no sensible, time dependent, definition of particle exists which will obey this 
criterion. This is related to the fact that particle production is stochastic and what is usually computed is the mean 
number of particles, which is the mean of a stochastic process; if the variance is not small compared to the mean, one 
needs to review the entire philosophy. On the other hand, we expect the notion of particles to be reasonably well 
defined if the quantum state is 'quasi-classical'. This brings us to the question of defining the notion of classicality of 
the state in some precise sense and relating it to particle production. 

The standard approach to understanding the evolution of the C—q system involves starting with the path integral for 
the full system and integrating over the q degree of freedom to obtain an effective Lagrangian L^ff = Lo(C) + Lcorr{^) 
in terms of C, that determines the transition amplitude between the asymptotic in and out vacuum states of the 
quantum subsystem when they are definable. This effective Lagrangian, in general, can have an imaginary part, 
which can be related to the mean number of g-particles produced out of the vacuum by the C-field over the entire 
evolutionary history of the system. (Roughly speaking it specifies the out-particle content of the in-vacuum state.) 
In order to compute the backreaction of the quantum subsystem on the classical degree of freedom, one needs to vary 
Lgff{C). However, since Le// can be complex, one would end up with complex equations (and solutions), which are 
difficult to interpret (see for e.g. [4]). What is usually done is to consider only the real part of the effective action 
in the calculation of backreaction and work with Lo{C) + KeLcorr{C). This term KeLcorr{C) is normally associated 
with vacuum polarization, and gives a contribution to the backreaction even in the absence of particle creation. On 
the other hand the vacuum-persistence-probability is directly related to the total number of particles produced and 
is determined by ImLcorr(C')- In such an approach it is not clear whether the back reaction due to the production 
of particles is incorporated in KeLcorriC). The dropping of the imaginary part, which also carries information about 
the quantum state of q, needs to be justified and one would like to have a formalism which provides a unified picture 
of the evolution of the quantum system in the background of C. 

We will now briefly connect up the above comments, presented abstractly in terms of a C and q, in two specific 
contexts widely studied in the literature. The first one is pair production in a constant electric field and the second 
one quantum field theory in an inflationary universe. 

In the case of the so called Schwinger effect in QED [5], we study the interaction of a quantized, charged scalar 
field (j) with a constant background classical electric field E (see refs. [6, 7] for a small sample of recent work and 
reviews) . Usually one computes the in-out transition amplitude by integrating over the ^ variable in the path integral 
[say, using Schwinger's proper time formalism] to obtain the one-loop effective action (known as the Euler-Heisenberg 
effective action) for E: 

(0+|0_)e= [ Vcj) e' J' ^^'^''-^^^^'^ = e'^^ff = e'/^-/^''*'". (1) 



The effective Lagrangian £e// turns out to be complex, and its (renormalized) real and imaginary parts are given by 



^ , ds 2 
ReCeff = - ^^^y I — cosm s 
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The real part possesses a Taylor expansion in and the lowest order correction to the classical Maxwell Lagrangian 
can be easily computed. The standard approach to computing the corrections to Maxwell's equations proceeds using 

only this real part; since it is analytic in one can obtain a modified action incorporating quantum corrections order- 
by-order. On the other hand, the density of pairs produced by the electric field can be obtained directly from the 
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imaginary part of the effective Lagrangian but is non-analytic in the couphng constant q. One would have expected 
the production of these particle pairs to drain the energy of the electric field and lead to a back reaction. It is not 
clear how an approach based on ReLef / encodes this effect contained in Re£e/ / and this needs to be clarified. 

As a second example, consider the standard inflationary cosmology. The quantity usually computed in this context 
is not the particle number density but is the magnitude of the power spectrum [8] . Since one uses it as a source in 
the perturbed Einstein equations, we are clearly moving into the domain of back-reaction. One normally computes 
the two-point correlation function of a quantized scalar field evolving in the inflationary background, and the power 
spectrum [the fourier transform in k-space] of the two-point function of the field is then used in calculating the 
fluctuations in density [9- 12]. This is different from computing an effective Lagrangian or from computing the number 
density of particles produced by the exponentially expanding universe (which is almost never attempted in the study 
of inflationary cosmology). The two approaches are thus technically different, and one needs to ascertain how the 
two-point function for the field is related to the particle content, and the extent to which the vacuum polarization 
part contributes to it. 

An added complication to this issue is that, although the early time evolution of the fluctuations [when the wave- 
length of a given perturbation is much smaller than the Hubble radius] is adiabatic and a choice of well-motivated 
initial conditions based on a natural definition of the vacuum state is possible, the late time phase [following Hubble 
exit, when the modes turn super-Hubble] is highly non-adiabatic and one can not identify a unique asymptotic out 
vacuum state. The notion of particles is ill-defined in such a dynamic background, and one needs to come up with 
appropriate variables that can help quantify the content of the quantum state of the perturbations in this background. 
Once again we need a generalization of our concepts from an out vacuum state to a time dependent situation. 

Another closely related issue stems from the fact that although the perturbations are assumed to be generated as 
vacuum fluctuations of a quantized field, their late time evolution at super-Hubble scales is treated using classical 
notions, with the only reference to their quantum origin being in the choice of initial conditions. Clearly, one needs 
to understand the validity of this assumption, and figure out how a quantum state can turn effectively classical. This 
has generally been regarded as an interpretational issue, and the usually adopted viewpoint in the literature hinges on 
the concept of decoherence [13] , which basically provides a mechanism for suppression of quantum interference terms 
between various states of the system (a generic quantum property) by taking into account its unavoidable interaction 
with a suitably defined environment. Mathematically, one proceeds by splitting the system under consideration into 
a part that is built out of that set of physical variables one is interested in, and an 'environment' that consists of the 
collection of those degrees of freedom that are in some sense inaccessible; the properties of the subsystem of interest 
can then be adequately described by a reduced density matrix that is obtained by starting with the full density 
matrix and tracing over the inaccessible degrees corresponding to the environment. Under suitable conditions, the 
reduced density matrix can be shown to reduce to a diagonal form (in a suitable basis) with the off-diagonal elements, 
representing quantum interference effects, getting suppressed. 

Although the formalism of decoherence can be invoked under fairly general conditions as real quantum systems 
seldom exist in isolation, ways have been suggested to do away with the quantum environment, and interpret classi- 
cality in terms of alternative notions. In the context of inflation, phase space correlations of the quantum variables 
have been analyzed using the Wigner distribution [10, 14, 15], and classicality interpreted in terms of peaking on the 
corresponding classical trajectory. The emergence of classicality has also been tied with the phenomenon of squeezing 
of the qTiantum state [10, 16]. Because of the potential implications of such approaches, which rely only on the 
intrinsic properties of the system for an understanding of quantum-to-classical transitions, there is a need to maJce 
them unambiguous. 

Keeping these questions in mind, we will conflne our attention to addressing two speciflc issues: (1) quantifying 
the physical content of a quantum state evolving in a time dependent background; in particular, identifying suitable 
variables that encode various aspects of the information contained in the quantum state, including the notion of 
particles, and (2) making the idea of interpreting classicality of the quantum state, as well as its relation with the 
concept of [appropriately defined] particle production, more precise. 

This will be done in a sequence of two papers. In the present paper, we will concern ourselves with studying the 
quantum mechanics of a single oscillator. This, of course, has general applicability. If one is studying a quantum field 
in an external time dependent background, like an FRW cosmological model [17] characterized by a scale factor a{t) or 
an electric field expressed in a time dependent gauge with vector potential A{t), the field can be decomposed by taking 
a spatial fourier transform into a set of uncoupled oscillators with time dependent parameters, each corresponding 
to a particular fourier mode k. This allows for a mode-by-mode analysis of the field evolution. Furthermore if, say, 
the oscillator labeled by the wave vector k is in the nth excited state, this is interpreted, in the field picture, as the 
presence of n particles with momenta k each. Presuming such a close correspondence with the field picture, we will 
freely make use of field theory language in our dealings with the quantized oscillator, with the understanding that the 
various definitions are to be looked at in the broader context of fields. Having said that, it may also be mentioned 
that when one considers quantum fields, some additional ideas can emerge. Firstly, the presence of an infinite number 



4 



of degrees of freedom brings in issues related to regularization, which need to be dealt with in a suitable manner. 
Secondly, one has the liberty of shifting from the oscillator picture by taking an inverse fourier transform back to 
coordinate space, and this opens up an alternative way to try understanding the physical content of the system. The 
case of quantum fields will be taken up in Paper II [18] . 

This paper is organized as follows. In section II, we outline a well-motivated and straightforward formalism, in 
the Schrodinger picture, to analyze the dynamics of a general quantized time-dependent oscillator, that proves a 
convenient tool to understand the content of the evolving quantum state, and the relationships between the various 
physically sensible quantities built out of it. Then, in section III, we first carry out a detailed analysis of a simple 
toy model, with particular emphasis on the question of relating particle content to classicality, and demonstrate that 
peaking of the Wigner function on the classical trajectory alone is not tantamount to having classical behavior; 
one needs to also look at some alternative measure of correlations to get rid of such ambiguities. We define such a 
reasonable alternative (calling it the 'classicality parameter') and based on an approximate asymptotic analysis of 
a class of models with general time dependent frequency Lu(t), show that whenever there is [appropriately defined] 
particle creation at late times, it is always accompanied by growth of the classicality parameter which can be quite 
unambiguously interpreted as an emergence of classicality. This is followed up with an analysis of a few additional toy 
examples, which serve to validate our analytic approximations and further illustrate the general features associated 
with the time evolution. Finally, we conclude in section IV with a discussion. In what follows, we shall set ?i = c = 1. 



II. EVOLUTION OF THE QUANTUM STATE: GENERAL FORMALISM 

The starting point of our analysis is a general one-dimensional harmonic oscillator with the following action, where 
we treat the mass and the frequency of the oscillator to be time dependent: 

A[q] = ljm{t)[f-u;\t)q']dt (4) 

with the conditions that m(t) > 0,0)^(4) > 0. (If m(t) is also monotonic, then using a time coordinate t with 
dr = dt/m{t) will get rid of m{t) and change to co'^m'^. But we will see that it is convenient not to do this and 
hence we will retain this form.) The classical trajectory of the oscillator can be obtained by solving the following 
equation of motion: 

j^(m{t)^^+m{t)ivHt)q = 0. (5) 

As mentioned earlier, the action in Eq. (4) can describe a particular fourier mode of a quantum field in the presence of 
a time-dependent classical background. In the standard procedure of canonical quantization, the variable q describing 
the coordinate of the oscillator is replaced by the corresponding hermitian operator q satisfying the commutation 
relation [q,p\ = i. 

In the Schrodinger picture which we shall work in, the evolution of the quantized oscillator can be described by a 
wave function ■ip{q, t) satisfying the time-dependent Schrodinger equation: 

.dtp{q,t) 1 d'^tl){q,t) 1 /,x 2/,^ 2 ta\ 

^^^ = -^^^+2™(^V {t)l ^M. (6) 

We are interested in the solution to this equation which could have been interpreted as the ground state of the 

oscillator at some time t = to. This suggests looking for a particular class of solutions of Eq. (6) that has been studied 
in various contexts [11, 14, 19 22], the form invariant Gaussian state, which is an exponential of a quadratic function 
of q. If one restricts oneself to states which have a vanishing mean (which happen to be the relevant ones in the 
examples we will consider), then the wave function assumes the form 

i;{q,t)=N{t)exp[-R{t)q^] (7) 

where N{t) and R{t) are complex quantities. Substituting the expression in Eq. (7) in the Schrodinger equation 
Eq. (6), the following equations for R and N are obtained: 

.N R . 2i?2 mu^ 

«T7 = — ; iR= —■ 8 
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Prom the first equation and its complex conjugate, it is easy to show that 



(9) 



which can also be obtained from the normalization condition on the wave function. Therefore, other than the overall 
phase of A'', the only non trivial aspect of the quantum state is encoded in the time dependence of the function R{t). 

Before we present the exact solution to this system, it is worthwhile to obtain the adiabatic limit of the system. If 
CO, m are slowly varying functions of time, then the Eq. (8) can be solved to give 



R{t) « m{t)uj{t) /2; N{t) = Nq exp 



(10) 



This determines the evolution of our state in the adiabatic approximation. Since R{t) « m{t)uj{t)/2 in the adiabatic 
limit, it is convenient to introduce a new complex function z{t) in place of R{t) by the relation 



R{t) 



muj fl — z 



(11) 



Clearly, z measures the deviation of R from the adiabatic value and we will call it the excitation parameter. The wave 
function is completely determined by z{t) as 



= Nexpl-Rq^] = AT exp 



muj ( 1 



1 + z 



A^exp 



mw ( 1 



lilmz \ 



1 



2Rez 



where 



|7V| 



2 mo; (1-^2) 



|l + ^r 



(12) 



(13) 



From the equation satisfied by R in Eq. (8) one can obtain the equation satisfied by z\ one finds that z satisfies a 
rather simple first order differential equation; 



z + li^oz - 



1 



TO 
TO 



(z2-l)=0. 



(14) 



Once this equation is solved for z with some appropriate initial condition (z = set at some instant to) we have 
completely solved the problem. Our analysis will be based on the time evolution of z. 

The evolution equation for z can be written in a slightly different form by using instead of the time coordinate 
t as the independent variable; given any monotonic range of ^(i), one can in principle invert this relation to obtain 
time as a function of the frequency, t = t{u)). With this replacement, Eq. (14) can be recast in the following form: 



dz(u}) 
to — — + 
duj 



^.(.) + l(.V)-i)=o 



(15) 



with e = {{d)/u)) + {rn/m))/ui, which will henceforth be called the adiabaticity parameter (and which too is expressible 
as a function of uj). This equation determines the function z{ui), and this conversion will be of particular relevance in 
our discussion of the effective Lagrangian in section II B. 

It may be mentioned in passing that the adiabaticity parameter as defined above has a simple physical interpretation 
in the context of quantum fields in a cosmological background characterized by the scale factor a{t). For a massless 
minimally coupled scalar field, any given fourier mode (describable by an oscillator) has a time dependent mass a^{t) 
and frequency [k|/a(/:), so the adiabaticity parameter is ek.(t) = 2a(t)/[k[, which (apart from a numerical factor) is 
just the ratio of the physical wavelength of the mode, Xp = 27ra|k|~^, to the Hubble radius Rh = {d/a)~^. 

The equations for R, z are first order but nonlinear. Since they are of the generalized Riccati type they can be 
transformed into second order linear equations. For example, if we set i? = — (zm/2) (/i//i) where fi{t) is a new 
function then Eq. (8) implies that /z satisfies the following differential equation: 



/i H fi + u) H 

TO 







(16) 
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which is same as the classical equation of motion, Eq. (5), satisfied by the oscillator variable q. The solutions of 
Eq. (16) obey the Wronskian condition 

W 

Hfi* - = (17) 
m(t) 

where W is independent of time. The variable z can now be expressed as 

<t) = 7^ • (18) 




As mentioned in the introduction, we now have completely solved the problem. One can, in principle, solve Eq. (16) 
to obtain the two linearly independent solutions for /i, say s and s* . The general solution is a linear superposition 
of the form /x = [As + B s*], with A and B determined by one's choice of initial conditions. Once has been found, 
the function R can then be computed, which completely fixes the quantum state of the system. The wave function 
tp{q,t) depends only on the ratio TZ = B/A, since it is independent of the overall scaling of /i. The real difficulty is 
not in doing this but in understanding what the physical content of the quantum state is. We shall now address this 
question. 



A. Particle content of the quantum state 

The quantum state given by Eq. (7) is more generally known in the literature as a squeezed quantum state and 
has been quite extensively studied [23] especially in the context of quantum optics. In the squeezed state formalism 
in the Heisenberg picture, one usually introduces the "squeeze" operator, defined as S{^) = exp [(l/2)(^*a^ — ^a^^)] 
with ^ = re*^, and where r is known as the squeeze parameter (0 < r < oo). The function z{t) that we defined earlier 
happens to be indirectly related to the squeeze parameters, as 2; = — e*^~^'''tanhr (where p = ui). Unfortunately, the 
squeezed state formalism does not help much in our tasks in hand; so we shall follow a different approach here. 

A physically motivated and reasonable way of quantifying the time-dependent content of the state would be to 
compare it with the instantaneous energy eigenstates at any given moment. Since the oscillator parameters are 
tini(^-dependent, one can not define stationary states in the usual manner; the alternative would be to define a set of 
instantaneous eigenstates at every instant. A wave function that starts off in the instantaneous ground state might, 
at a later time, be in a superposition of instantaneous eigenstates defined at that moment; this can be thought of as 
excitation of quanta. 

We will define these instantaneous eigenstates, at a given moment t, as a. set of states that have been obtained by 
adiahatically evolving the eigenstates at some initial instant to. This amounts to arranging matters in such a way that 
if the oscillator evolves adiahatically, a wave function that has been set in the instantaneous ground state a,t t = to 
will evolve to coincide, at every subsequent moment t, with the instantaneous ground state defined at t. This set of 
instantaneous eigenstates at time t has the form 



, , , /mw\ V i 

<Pn(q, t)=[ ^ — HnWmujq) exp 



mu o . /"* / 1 



to 



2 q^-i j {n+^-]uj{t)dt 



(19) 



where Hn are the Hermite polynomials [24] and n = 0,1,2,.... Except for the phase factor, these are essentially the 
harmonic oscillator wave functions, constructed using the instantaneous values of Lo(t)^m{t) at time t. They form a 
complete set of orthonornial functions at time t and any other wave func;tion can be expanded in terms of them. The 
phase in Eq. (19) is so chosen that it coincides with the phase of the wave function in the adiabatic limit obtained in 
Eq. (10). 

We can now expand our exact solution -f/; in terms of the instantaneous eigenfunctions. Since if) is an even function, 
the amplitude for the oscillator to be in the nth instantaneous eigenstate (/>„(t) at time t is non-zero only for even n, 
and is given by 



Cn{t) = r dq4>l{q,t)i,{q,t) = N (^)'^' T dgif„(V^9)e-(^+^)'^'+'^4("+^)-^^^-. (20) 

The above integral can be evaluated by standard methods (see Appendix A for details). We can then compute the 
probability associated with the amplitude C„(f) to obtain 

i'2n[t) - \C2n\ - ^J^^^ (21) 



{t)dt 
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with the time- dependent vacuum-persistence-probability Po(i) given by: 

Po{t) = P= ^/l^W (22) 

Clearly the probability for occupying the excited states is controlled by z justifying the name excitation parameter. 
The generating function for the probability distribution in Eq. (21), defined as G(x) = Yl'^=o ■^^nX"' , is given by 

Once the probability distribution is known, it is trivial to compute the mean number of quanta in the state at any 
time t; this is given by 



in) = ^ 2nP2„ = 2 G'(l) = tj^^ = (24) 



We shall take this quantity, (n), as describing the time dependent 'particle' content of the quantum state. We do 
not claim that our definition is unique in any sense; only that it will be useful and physically reasonable. This 
interpretation will be borne out by different considerations in what follows and the first of these is the following. The 
system evolves under the action of a time-dependent Hamiltonian H and one can compute the mean value of the 
energy at any time t by the expectation value of the Hamiltonian, E{t) = {tplHltjj). Direct computation shows that 



= (2^) (1-^1' + '^' I'^H ( + 2 ) ^^^^ 



This clearly strengthens the motivation to think of (n) as defined above as the number of quanta present at time t. 

Let us rewrite the expression for the vacuum persistence probability, Eq. (22) in terms of the the mean particle 
number: 



Po{t) = vT^N^ = (1 + («))"'^' = exp 



+ in)) 



(26) 



It is obvious from this expression (as well as, of course, from Eq. (21)) the excitation probabihty for different levels 
is not Poissonian. When the excitation to level 2m is interpreted as creation of m pairs of particles in quantum field 
theory, this implies that pair production is not — in general — a Poisson process and is non-trivially correlated. In 
the limit of (n) <C 1, however, we have Pq ~ exp(— (n)/2) which allows one to identify (n)/2 as the mean number of 
pairs produced from the vacuum, which is a standard result. 

It is clear from the expression for the generating function, Eq. (23), that the mean particle number (n), as well as 
the higher order moments, are functions of only the magnitude of z. On the other hand, the wave function in (12) 
is built out of not only the magnitude, but also the phase of z. In fact, writing z = \z\e^^ and using Eq. (24) to 
determine \z\, we find that the wave function can be expressed in the form: 



V'(*>9)=(— j (^1 + 2(n) + 2^{n){{n) + 1) cos 6lj exp 



1 2/ l-2iJ(n)((n) + l)sm.9 
■-mwq = 

2 \l + 2{n) + 2^{n){{n) +1) cos 9 



(27) 

It is obvious that the particle number only contains partial information about the quantum state; given the particle 
number (and even all the higher moments) at a given moment of time, one can not have complete knowledge of the 
state of the oscillator at that moment. 

While on this issue, it is worth noting that (n) 00 when \z\ 1. In this limit the width of the wave function 
scales as (n) and the gaussian spreads all over the coordinate space. In the same limit the width of the wave function 
in momentum space goes to zero, "squeezing" the wave function to the p = axis. One can also show that the 
dynamical equations are same as the following equations for (n) and 9: 

^n)=(^+'i^) v/(n)((n) + l)cos^; 9 = -2co - + -M±^sin9. (28) 



It is obvious from these equations that even when lu and m have a fixed sign, the sign of (n) depends on the phase 
9 and hence need not be monotonic in general. Moreover, it may be noted that if (n) is specified as a function of 
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time, one can, using the above equations, determine the form of 0(t), and this allows one to reconstruct the complete 
time evolution of the wave function. In this sense, it is possible to fix the state of the system completely given the 
time-variation of the mean particle number alone. And further, if (n) is a monotonic function of time, then one can 
trade off the t dependence for dependence on (n), and thus express z (as well as the wave function) explicitly in terms 
of just the mean number of particles. Such a transformation will be well-defined only when there is a one-to-one 
relation between t and (n) . (If the mean particle number is oscillatory, then there would in general be several values 
of time for the same (n), and consequently the expression for the wave function in terms of (n) will be multiple- valued 
as well, and not uniquely specifiable.) 

Getting back to our quantum state, it is clear that an additional physical variable needs to be specified that can 
provide information about the phase 9. A suitable choice for this is the spread in the wave function, given by the 
expectation value of q'^. It has the form 

2 



(,^)=£,"W<,.t)l"*=5^ = |L. (29) 



In the context of quantum fields in a cosmological setting, this dispersion can be directly related to the logarithmic 
power spectrum [9, 10, 12, 14, 19], which is the fourier transform of the two-point correlation function for the field 
evaluated in a particular state: 

'''Pt'ik) = ^{ql). (30) 
The quantity (q^) also can be re-expressed in terms of 2; as follows: 



2mu 



z 



2 



(2(n) + l) , ((n) + l) 



+ ' ^ Re{z). (31) 



2moj mu 



This expression, again, shows that the power in the q-mode is not completely expressible in terms of the instantaneous 
mean particle number [which involves just the magnitude of z], because it carries additional information encoded in 
the phase of z. It follows that, in general the computation of the power spectrum in Eq. (30) (in the cosmological 
context, for example) is different from computing the mean number of particles produced. However, it can be seen 
from Eq. (31) that in the limit when the real part of z becomes a constant, {q^) can be written purely as a function 
of the mean particle number. 

Our treatment is completely ecjuivalent to the standard analysis that is generally done in the Heisenberg picture, 
in which the time dependence of the system is encoded in the Bogolyubov coefficients a{t) and /3(t) satisfying the 
relation |a(t)p — |/3(t)P = 1. The function z that we have defined is related in a simple manner to these Bogolyubov 
coefficients, as 

z{t) = ^ e-^'P^'^ (32) 

with p{t) = u){t) (see Appendix C). In fact it can be shown quite trivially, that if one sets the oscillator in the vacuum 
state at time to, then the above quantity is directly proportional to the [normalized] amplitude for the vacuum state 
defined at time to to be a 2-particle state with respect to the vacuum defined at time t: 



{2,t I 0,to) ^ 1 (0,t \a{t)a{t)\ 0,to) ^ J_0Htl ^ ^ 
{O,t\0,to) V2 {O,t\0,to) V2a*{t)^^' 



2ip{t) 



(33) 



where a{t) is the annihilation operator defined at time t, with a(t)|0,t) = 0. It follows from the relation in Eq. (32) 
that all the physical quantities are alternatively expressible in terms of the Bogolyubov coefficients. In particular 
Eq. (24) now reduces to (n) = which is the standard result. We however shall not be using them in the analysis 
that follows, but will stick to understanding the evolution of z{t). (In field theory, each fourier mode labeled by a 
wave vector k will have a corresponding Zk and one can also obtain the spatial fourier transform of this quantity; we 
will see in Paper II [18] that it contains valuable information about the classicality of the state.) 



B. Effective Lagrangian 



Another useful way to quantify the effects induced by the time-dependent background on the quantum system, 
especially in the semiclassical context, is to look at the vacuum persistence amplitude, which measures the amplitude 
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for a state to be a vacuum at late times, if it started out as a vacuum state at early times. Normally, this quantity is 
defined using asymptotic in and out vacuum states [1, 2, 5]. It is directly related to the effective action, the imaginary 
part of which specifies the asymptotic particle content in the quantum state, while the real part is used in analyzing 
backreaction. 

This idea needs generalization when asymptotically adiabatic vacuum states cannot be defined. One would like 
to have some sort of a generalized time-dependent analogue of the effective action that is amenable to a suitable 
interpretation. Based on the formalism outlined in the previous section, one can define such a quantity in a fairly 
natural manner. 

Consider a situation in which the oscillator has been set in the instantaneous ground state at some instant, say, 
t = to. This fixes the form of z{t) and hence of R{t). One can then compute the amplitude for this oscillator to 
evolve and be in the instantaneous ground state at some moment t > to in the future: this is just the amplitude Co 
evaluated in Eq. (20), and one can write it in terms of a time dependent 'effective action' as follows: 



N{l) [iii.'^'^^^'^ 



Co{t) = " y ^ V" — "/ ^ j ——(jn = exp{iAe//(f)} = exp i / Leff{t)dt 

Jtn 



to 2 



^{Rit) + ^) 

Since the amplitude is known it can be easily computed and turns out to be a rather simple expression: 



(34) 



i^e//(i)=^ + (35) 

To place this result in context, recall that the time dependence of ui, for example, arises because the classical degree of 

freedom (background metric, electric field ) is described by a time dependent solution. This makes quantities like 

w = {duj/dC)C explicit functionals of C. Further, the effective Lagrangian also depends on C{t) implicitly through 
z which is determined in terms of the background variables through the differential equation Eq. (15). While the 
effective Lagrangian looks rather simple when expressed in the above form this simplicity is deceptive (and useless 
for the purpose of calculating backreaction etc.). The effective Lagrangian needs to the thought of as a functional 
of background variables and this is often nontrivial. Nevertheless, we can make significant progress using the above 
expressions, keeping the relevant caveats in mind. 

The real and imaginary parts of the effective Lagrangian are directly related to the function z. The real part can 
be written as 

ReLeffit) = -^(^^ + ^^lmz=-^eu) Imz. (36) 

Note that we have defined the particle content and the vacuum state using the instantaneous eigenstates obtained by 

adiabatic evolution. Therefore, the usual adiabatic term (integral of oj over time) does not appear in this expression 
and ReLgff vanishes in the adiabatic limit (e ^ 0). The imaginary part, on the other hand, is given by 

lmL,jf{t) = U- + -)Rez. (37) 

We would have expected it to be related to the particle content and indeed it is. Using the equation for (n), one can 
show that this is equal to 

ImLe//(t) = ln(l + (n)); ImAe// = ^ ln(l + (n)) (38) 

This is in accordance with the standard interpretation of A^ff if the system has a late time adiabatic regime, where 
lmAeff{t) saturates to a constant value specifying the asymptotic particle number. If (n) <C 1 then the imaginary 
part of the effective action becomes ImAeff{t) ~ (1/4) (n) so that the vacuum persistence probability is 

|(0,i|0,io)|' « exp(-2ImAe//(t)) « 1 - ^(n) (39) 

which matches with the result in Eq. (22) for |zp <C 1. (The factor (1/2) is due the fact that (n) is the mean number 
of particles while in quantum field theory one usually quotes the result for mean number of pa/i/rs.) We will see that 
our definitions make sense in all other cases too when we look at the behavior of L^f f in different limits in the study 
of toy models in the next section. 
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To recapitulate, we have defined the quantum state in terms of a parameter z and find that the particle content of 
the state determines (and is determined by) . To fix the state uniquely we also need to know the phase of which 
can be obtained from some other suitably defined quantity like the dispersion l^q^). Together they completely specify 
the state of the system at any moment. The next question we want to address is the 'classicality' of this state, and 
its approach to classicality. This may be understood, in one possible way, by shifting attention to the system's phase 
space. 



C. Wigner function 

We now will attempt to quantify the level of classicality of the quantum state in terms of its phase space correlations. 
A suitable tool for this purpose is the Wigner distribution function [14, 15, 19], which is defined, for a wave function 
V'(g,i), as 

OO 

= ^ I dur{^^\i) e'""- (40) 

— OO 

The Wigner function can be regarded as a quantum analogue of the classical distribution function, and satisfies the 
following identities: 

/•oo 

W((z,p)dp= |V(g)P, (41) 



J — c 



y^{q,v)dq=\^{v)? (42) 



where ^piji) is the wave function in Fourier space (i.e, the fourier transform of V')- These two relations suggest that 
yV'((j',p) can be thought of as a probability distribution in phase space, provided it is positive. In general, Wigner 
function satisfies the evolution equation: 

5W paw dV dYV h"^ d^V d^YO 
dt m dq dq dp 24 dq^ dp^ 

where V{(i) denotes the potential, and dots indicate terms with higher powers of h and higher derivatives of W and 
V. For quadratic potentials, the right hand side of eq.(43) vanishes and we recover the classical continuity equation. 
(This is true for any potential up to order fi^.) Although the above relations suggest that the Wigner function might 
be interpreted as a 'joint probability distribution', it can take on negative values in some cases [23]. Gaussian states, 
however, turn out to have a positive-definite Wigner function. 

The Wigner function is useful in studying quantum to classical transitions in a system in terms of correlations 
between the phase space variables {q^p); a pure quantum system is represented by a completely uncorrelated Wigner 
function like W{q^p) = A(q)B(p) so that the probability for the system to have a momentum p is independent of its 
position q. The ground state of the harmonic oscillator, for example, has such a Wigner function showing it is very 
non-classical. For a classical system, we will find W to be peaked in a limited region in phase space with the totally 
classical state being W{q,p) cx S{p — ,f{q)), where p = f{q) represents the classical phase-space trajectory. Hence the 
classical system is expected to show a strong correlation between q and p. A useful measure of the 'classicality' of a 
state can thus be provided by the degree of correlation between q and p. For this purpose, we consider the following 
quantity: 



with the average for any function F{q, p) of the phase space variables being calculated using the Wigner function: 



oo /.oo 



{F{q,p))w^ / F{q,p)W{q,p)dqdp. (45) 



oo J —CO 



We shall henceforth refer to the object S as the classicality parameter, for a pure quantum state, like the ground sate 
of the oscillator, we get 5 = 0. For a highly classical state, we will expect S ^ 1. 
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The Wigner function corresponding to the gaussian wave function in Eq. (7) can be expressed as 



yV{q,P,t) = - exp 



a'{t) {p-J{t)qY 



where a and J are given by 

= {R + R*)~^ and J = i{R-R*). 
The above expressions can be written in terms of the function z (and (n)) as follows: 



^2 ^ 2|/z 



\l + z\ 



and 



J 



W mw(l — I^P) mto 
m d(Zn|/tip) 2mujlm{z) 



2 dt 



For this Wigner distribution, the classicality parameter is expressible in terms of cr^ and J: 

^/l + (Ja2)2■ 



(46) 



(47) 



(48) 



(49) 



(50) 



It can be seen from eq.(46) that when J = the Wigner function represents an uncorrelated product of gaussians in 
q and p, and 5 = 0; this happens to occur when the gaussian state coincides with an instantaneous ground state. On 
the other hand, when both cr^ and J take on non-zero values, W(g,p) becomes correlated and iS would be appreciably 
different from zero. [It may be noted that |»S| < 1, and the maximum possible value that it can have is unity.] The 
choice of 5 as a measure of classicality is thus a reasonable one. 

The idea of associating particle creation with the approach to classicality of a quantum state has been around in the 
literature for quite some time, particularly in the context of cosmology [14, 16] . It has been demonstrated, in the case 
of perturbations evolving during an inflationary epoch in the early universe, that the Wigner function gets peaked 
on the corresponding classical trajectory as the state undergoes large squeezing at super- Hubble scales [14, 19]. (The 
argument in this case hinges on the fact that the frequency of the oscillator turns imaginary at late times, and results 
on peaking of the Wigner function for the inverted oscillator are then applied to draw the relevant conclusions.) Here, 
we are interested in exploring the generality of this connection. In particular, we would like to know how production 
of partic;les (identified, in our c;ase, with a growth in the quantity (n)) is related to the spreading of the Wigner 
function, as well as to the variation in the sharpness of the q-p correlation, i.e. the classicality parameter. 

We note that for the time dependent oscillator (with completely general m{t) and a;(t)), both the particle number 
and the q-p correlation can be built out of the complex quantity z; but there is no simple relation directly connecting 
the two. The classicality parameter involves phase information of z as well, and this makes drawing conclusions 
about its behavior solely on the basis of one's ideas about (n) (which depends only on the magnitude of z) not 
so straightforward. This means that one needs to actually determine z to address the question of relating the two 
variables, to which we now turn. 



III. ANALYTIC APPROXIMATIONS AND ASYMPTOTIC ANALYSIS 

We move on to explicit analysis of several concrete examples based on the ideas of section II. For convenience, 
we set m = 1, so that the time dependence enters only through uj{t). (As we explained earlier, this entails no loss 
of generality.) The appropriate measure to characterize the behavior of the oscillator is, of course, the adiabaticity 
parameter given by e{t) = lo/lo"^, whose magnitude determines the nature of its dynamics. We will assume the late 
time limit to be given by the time variable t going to infinity, and consider the two extreme cases: one in which the 
late time evolution is adiabatic (e <C 1) and the other in which the adiabaticity condition is strongly violated (e 1). 

For the case of an oscillator with a time-dependent frequency oj{t)^ the equation (14) for the function z, simplifies 

to 

z + 2iu)z + ^{z^ - I) =0. (51) 

As mentioned before, once z is known, all other variables of interest can be trivially obtained. The table given below 
encapsulates the expressions for the various physical quantities in terms of z for ready reference: 
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The wave function 


^-(9,*) =exp[-ijg2] 


Evolution equation for z 


z + 2iuz + ^ (z^ - 1) = 


Mean particle number 




Wigner function W {q, p, t) 
= i exp[-^-<T^(t) ip-J{t)qf] 


^ _ 2u,/m(») 


Classicality parameter 
5 — ^"^ 


2Jm(») 


Spread in the wave function 
(related to the power spectrum) 




Effective Lagrangian 


ReL^ff — — ^ ^ Imz 
ImL^jf = L!kRez=i-^ ln(l + (n» 



We axe particularly interested in how an oscillator, that starts off in the instantaneous ground state at some moment 
t = to, would evolve at late times [the t ^ oo limit]. With this in mind, wc will specifically consider a scenario in 
which the system enters an adiabatic phase at late times and, alternatively, when the asymptotic evolution deviates 
from adiabaticity. All the useful choices for the frequency ui{t) fall within one of these categories. Before we move on 
to carry out an approximate analysis for such functions, we will take up a simple toy model first. 



A. The case of constant adiabaticity parameter 

We will begin our analysis by considering the case of constant adiabaticity parameter e. (In the cosmological 

context, this model would describe a massless scalar field evolving in a background characterized by the scale factor 
a(t) oc t, which corresponds to having a matter source with equation of state p = —p/'-i. This is the borderline case 
between accelerating and decelerating models.) This model is exactly solvable, and can describe adiabatic as well 
as non-adiabatic evolution through appropriate choices for the value of e. This flexibility allows one to see how the 
nature of the evolution changes with transition from adiabaticity to non-adiabaticity. This example is also expected 
to shed some light on the general features one may expect to find in the two extreme limits. 
For the case of constant e, the frequency function is given by 

.[t) = (0 < i < oo) (52) 

with J/ > 0, A having the dimension of frequency and the adiabaticity parameter being equal to —v. The combination 
(j/A)~^ sets the time scale for variation; for t » {v\)~^, the frequency falls as l/{vt). For this w(f), equation (51) for 
z can be solved analytically to obtain the following result (where t = Xt and we have set 2; = at t = 0) : 

Zi[l + VT) 2 - Z2 

with zi,2 = (2? ± Vi^^ - 4)/z/. 

It is convenient to analyze the behavior of this model by splitting the range of u into two parts: one corresponding 
to z^^ > 4 and the other to v"^ < 4. Let us consider the former case first. For > 4, the expression for the particle 
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(a) (b) (c) 



FIG. 1: Variation of the mean particle number (n), classicality parameter 5 and the excitation parameter z with r for constant 
|e| = 20. (n) is a monotonically increasing function of time. S increases sharply and saturates at unity at large times. The 
trajectory of z remains confined within a quadrant and ends up at a point with \z\ — 1 in the r — > oo limit. 



number is given by 



This is an exact expression valid at all times. It is clear that the particle number increases monotonically with time 
in this case. Since we are particularly interested in the scenario of highly non-adiabatic evolution, it is appropriate to 
consider large values of v (in comparison with unity). For i/ 3> 1, we have the following approximate expression for 
the particle number: 

1 I^^T^ 

Let us consider the early and late time limits of the above expression. For the case of i^r <^ 1, 

(n)«^^ + 0(^V). (56) 



On the other hand, for i/t ^ 1, 



VT \ , ^ 



For strongly non-adiabatic evolution, starting from a vacuum state, the particle number increases without bound at 
late times. But Eq. (57) also shows that the energy content of the particles (n)a;, however, saturates at a constant 
value (1/4) at late times. Our approximation is borne out by the plots of the particle number and the quantity {n)Lu, 
shown in figures 1(a) and 2. The particle number grows without bound, while the latter quantity has a constant 
limiting value. 

The classicality parameter S for this case, describing the strength of phase space correlations, is plotted in figure 1(b). 
It is evident that S starts from zero [corresponding to the initial vacuum state], and quickly grows to unity as 
the particle number increases with the progress of time. This clearly presents an example where the phase space 
correlations grow in accompaniment to particle creation due to violation of adiabaticity. 

One can also directly visualize the trajectory of the state in the complex z plane, which has been plotted in 
figure 1(c). z starts from the origin, corresponding to the initial vacuum state, stays within a quadrant, and ends up 
at a limiting point corresponding to |z| = 1 at late times. 

One would also like to know how the Wigner function evolves with time. This can be understood by computing 
the functions and J using the relations (48) and (49). These functions have been plotted in figure 3. The value of 
starts from l/cij(i — 0) and monotonically increases with time, while J starts from zero, rises to a maximum, and 
then falls off again to zero at late times. This variation corresponds to a Wigner function that ends up peaking on 
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FIG. 2: Variation of {n)u! with time for the case of |e| — 20. This quantity saturates asymptotically, so the energy remains 
finite at late times. 





FIG. 3: Variation of the functions and J for |e| — 20. grows monotonically while J' falls to zero at late times, indicative 
of a Wigner function that peaks on the axis. 



the g— axis at late times. Let us compare this with the way the corresponding classical trajectory behaves with time. 
The classical trajectory is given by 



qdr) = Vl + 



VT \ cie 



and 



Pc{T) 



f^) (1 + + ( 1 - f7I\ (1 + 



(58) 



(59) 



where ci, C2 are real constants. It is clear that in the late time limit, qc oo while the momentum ^ 0, so every 
trajectory at sufficiently late times ends up on the q—axis, coinciding with the peaking of the Wigner function. This 
feature, together with the behavior of the classicality parameter, provides strong indication of a quantum-to-classical 
transition taking place at late times. 

The numerical value of the effective Lagrangian for this case can be evaluated using Eq. (35). (As we mentioned 
before, such an expression has a limited validity since the effective action should be treated as a functional of w rather 
than as a function of t for a particular LLj{t). Nevertheless it illustrates certain interesting features.) In the highly 
non-adiabatic limit of 3> 1, it has the approximate form 



A 



2(i/T + 2)2 4 iiyT + 2) 



This reduces to the following limiting form at early times [corresponding to vt 1]: 



(60) 



(61) 



Interestingly, the real part of the effective Lagrangian can be re-expressed in terms of the particle number using 
Eq. (56): 



ReL, 



{n)uj/2 



(62) 
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FIG. 4: Variation of the mean particle number, classicality parameter and z for constant |e| = = 0.2. Both(n) and S are 
oscillatory at all times with constant finite amplitudes. The complex trajectory of « is a circle in the upper half-plane, reflected 
in the oscillatory nature of (n) and <S. 



This shows that the real part of the effective Lagrangian does contain information about the particle production (and 
hence can encode it in the back reaction calculation) in the non-adiabatic, but (n) -C 1 limit. We will find later that 
this is indeed a general feature. On the other hand, in the late time limit, 



A 



1 iv 1 

IVT 4 VT 



A 



1 i V 
(n) + 2(n) 



for VT » 1. 



(63) 



Prom the above expression, it is clear that although Lpj^ij) ^ in the late time limit, the effective action is 
logarithmically divergent, since ^e// oc t~^. This is related to particle production proceeding without bound at late 
times. 

We next consider the other regime, covering values of the adiabaticity parameter lying in the range < 4. For 
this case, Zi , are purely imaginary and the expression for z may be rewritten as 



iy/^-llnCl+i/r) 



z\e 



y^-lln(l+i/T) 



The particle number has the exact form 

(n) = 



(4 - z/2) 



sm 



lln(l + z/T) 



(64) 



(65) 



This quantity keeps oscillating between the values and ji^— i^^), representing a continuous interplay of creation 
and annihilation of quanta. This situation thus presents a case where the oscillations persist with constant amplitude 
and the particle number never settles down to a constant value [or enters a phase of monotonic variation]. We have 
plotted the variation of the particle number as well as the classicality parameter for this case [with v chosen to be 
0.2] in figures 4(a) and 4(b). The classicality parameter is also oscillatory in the absence of a monotonic variation in 
(n), and this may be contrasted with the non-adiabatic case considered earlier, where iS grows and saturates at unity. 
This is suggestive of the possible close connection between particle creation and growth of phase space correlations, 
which will be explored in greater depth in the examples that will follow. 

The complex trajectory of z(t) is also depicted in figure 4(c). The values of the real and complex parts of z 
oscillate with time, and the trajectory is a circle with the center displaced from the origin. The magnitude of 2, too, 
oscillates and this is reflected in the variation of the particle number. This evolution differs sharply from that for 
large adiabaticity parameter, where the number of particles as well as the q-y correlation display monotonic increase 
with time. 

As for the Wigner function, its evolution is described by the functions and J ^ shown in figure 5. While J goes 
to zero at late times, is an increasing function of time, driving the Wigner function to get concentrated on the 
g— axis at late times. In comparison, the classical trajectory has the general form 



(66) 
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where Ci is real and </> is a constant phase factor, and 



Pc{r) = q'c{r) = 



2^1 + VT 



cos 




1 ln(l + vt) + 




1 sin 




- 1 ln(l + vt) + > 



(67) 



The ampUtude of the momentum steadily decreases, while that of Qc grows with time. The trajectory, consequently, 
ends up on the g— axis asymptotically. This matches with the late time behavior of the Wigner function. Thus 
the Wigner func;tion peaks on the classical trajectory at late times, similar to that found in the non-adiabatic case 
earlier, in spite of the fact that here there is no genuine particle creation. It may also be noted that neither does 
the correlation measure S grow monotonically at large times. This example thus brings out a potential problem with 
regarding peaking of the Wigner distribution on the classical trajectory as the sole indicator of classicality. At the same 
time, the behavior of the classicality parameter, which evidently tracks particle creation, and is markedly different 
in the two cases, suggests that this measure may be able to provide a useful means of resolving such ambiguities 
regarding the interpretation of classicality. 

In the strongly adiabatic limit, we have v <^1, and the particle number is approximately given by 



(n) 



sin^ -ln(l + z/r) +0(i/*). 



For small times (r 1), ln(l + vt) « vt and sinr w t, so that 

(n) 



9 • 2 

V sm T 



(68) 



(69) 



This behavior matches with the early time variation of (n) in the extreme non-adiabatic case (corresponding io v ^ 1). 
The effective Lagrangian for v^ < A takes the following form: 



ImLeff = - 



4 ((n) + 1) 



1 - 



sm 



-lln(l + z/T) 



l-cos(^^-lln(l + i.r)) 



(70) 



and 



ReLe// = 



LOV 



sm 



1 ln(l + vt)^ 



f (l+cos(y;^ln(l + i/r))) 



(71) 



In the adiabatic limit oi v <^ 1, the real part is approximately given by 



Rei 



e// 



utv'^ 4(n) 
~8 ^ 



cos 




1 ln(l + vt) 



uj{n) 



[1 + 0{v^ 



(72) 



and this is. of course, valid at all times. The identical form for KcL^jf in both the limits of small as well as large 
values of the adiabaticity parameter, which may be noted from Eqs.(61) and (72), suggests that when the particle 



17 



number is small, the real part of the effective Lagrangian is directly related to the particle content in the quantum 
state; this is of relevance to the issue of accounting for back-reaction due to particle creation using HeL^ff , a point 
we shall return to later. 

To summarize, the model of constant adiabaticity parameter, for different choices for the magnitude of e, can 
exemphfy both adiabatic and non-adiabatic evolution. In the adiabatic Umit of small |e|, particle creation is suppressed 
with the mean particle number and the classicality parameter (both starting from zero corresponding to an initial 
vacuum state) remaining finite and oscillatory at all times. This behavior sharply differs from that in the extreme non- 
adiabatic case of large |e|, where the mean number of particles as well as the q-p correlation diverge with time, although 
the energy saturates at a finite limiting value. In both cases, however, the Wigner function ends up peaking on the 
corresponding classical phase space trajectory at late times, in marked contrast to the behavior of the classicality 
parameter. 

Having gained some idea of the kind of features which may be expected in adiabatic as well as non-adiabatic 
evolution, we now move on to consider the case of a general adiabaticity parameter e{t), and examine its two possible 
late time limits. 



B. Adiabatic evolution at late times 



For a general uj{t) we can physically distinguish the following situations. Depending on the numerical value of 
e{t) one can describe the evolution as adiabatic or non-adiabatic. Further, in the c;ase of adiabatic evolution one c;an 
further distinguish between two cases, the first case corresponds to the system evolving in an arbitrary fashion till 
about t = T, say, and evolving adiabatically for t > T. In this case the quantum state will not an instantaneous 
vacuum state at t = T and certain amount of particle production would have taken place by then. The further 
adiabatic evolution should preserve this particle content with a higher order contribution to small amount of further 
particle production. The second case corresponds to an adiabatic evolution all the way from to to t with the system 
starting at the vacuum state at time t = t^. In this case, the particle production will be a higher order effect and we 
will be interested in the lowest nontrivial order contribution. We will first analyze the case of the evolution proceeding 
adiabatically at late times. This covers all possible choices for the function u){t) which fall off slower than 1/t in the 
f — > 00 limit, and correspondingly, e{t) 0. 

A perturbation theory approach can be adopted to find an approximate solution in an adiabatic regime. We start 
by writing the function f{t) in terms of two new functions A{t) and B{t) as 

/(f) = 4^e*e(*) + Me-'^(*) = u) (73) 

\/2a; \/2a; 



and imposing the additional constraint condition (to take into account our introduction of two new functions in place 
of the original one) that 

/(f) = ia;4^e'^W - ia;Me-»^(*). (74) 
This gives (details can be found in the Appendix C) 

— lUJ ' 



2ie(t) - 

(75) 



/ \A{t) + B{t)e-'^i'^(i) 
and yields the following equations for A, B: 

^=il^e-2ie« , ^=A^e2^^W. (76) 

The first order adiabatic approximation for the functions A{t),B{t) can be found by assuming initial conditions 
A{tQ) = a,B{to) — b (with a, b being arbitrary complex constants) at some initial time to within the adiabatic regime. 
To a first approximation, the coupled equations then give 

lb 



A{t) « a + z^(e(f)e-2^^(*)-e(io))-| / rfi'e(Oe"'*'^*'\ 

J to 

B{t) « b-i^ (e(t)e2'^« - e{to)) + '-^ f dt'e{t')e^'0^''l (77) 

J to 
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At the lowest order, one can ignore the integrals over e on the assumption of adiabaticity. (This result, therefore is 
exact for the the e = constant case discussed earlier.) So, at the first order of approximation, we have 



z{t) 
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A{t) 



-2ie(t) 



-2i0{t) 
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e(io)(l+^ 



-2ie{t) 



e{t) 1 + -oe- 



-4i0(t) 



+ 0(e2) 



(78) 



where the function 71(1) is of order e. This solution is valid in the period when the evolution is adiabatic, irrespective 
of the past history of the system. For example, one can envisage an evolution in which the system starts at the ground 
state, evolves very non-adiabatically for a period of time, ti < t < T say, and is adiabatic for t > T. In that case, the 
system will be in a highly excited state (with some amount of particle production already having taken place during 
ti < t < T) dX t = T and during t > T, this state will evolve adiabatically. In such a case we cannot say an}rthing 
about the magnitude of a,b and we will keep them arbitrary for the moment. (This situation should be contrasted 
with another one frequently discussed in the literature in which it is assumed that the evolution is adiabatic throughout 
the period. We will consider this case later on.) The corresponding \z\^, to order e, is given by 
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Using this expression in Eq. (24), one can compute an approximation for the mean particle number, which, again to 
order e, is given by 

n2 
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(80) 



Let us now consider various limits of this expression. 

If the system evolves very non-adiabatically for a period of time ti < t < T and is adiabatic for i > T then we 
can take the e(t) limit at t > T. Then, for any non-zero value of b/a, the terms dependent on e{t) in the above 
expression get progressively suppressed. This implies that the mean particle number has a finite limiting value given 
by 
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(81) 



This result is completely understandable. In this case, all the particle production takes place during ti < t < T and 
the constants A, B play the role of Bogolyubov coefficients for the evolution in Heisenberg picture; the mean number 
of particles is indeed If one compares the expression for z in Eq. (78) with Eq. (32) relating z to the functions 
a and /? in the Heisenberg picture, it can be seen that in the e{t) — > limit, the ratio of the Bogolyubov coefficients 

becomes a constant: (3* /a* w {b/a) + (ze(to)/4)(l + (&^/a^)). Using the Wronskian condition jap — = 1, the 
limiting expression for (n) derived above can then be shown to be equal to |/3p. This is exactly the expression one 
expects to obtain for the asymptotic particle number in an adiabatic 'out' vacuum state in the Heisenberg picture. 
This thus demonstrates the close correspondence between our approach and the conventional one. 

An expression for the q-p correlation can also be worked out, using the approximations in Eq. (78) and Eq. (79), 
to first order in e: 
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In the late time limit (as e{t) 0), the function TZ has a constant amplitude, and since 9 is an increasing function 
of time [owing to our choice for w(f)], the above expression is oscillatory with a finite amplitude. Thus, in the late 
time adiabatic regime, once the particle number has saturated, the phase space correlation {qp)w remains bounded 
and oscillates around zero, again indicative of the necessity of continuing particle excitation for driving the system 
towards classicality. 

One can also look at the effective Lagrangian, defined in Eq. (35), in this limit. The imaginary part of the action will 

become a constant at late times with its value related to the asymptotic number of particles produced, in conformity 
with the usual interpretation. The real part, using Eq. (36) and Eq. (78), is given by 
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to second order in the adiabaticity parameter. Written in this form, it is evident that the first 0(e) term is a rapidly 
oscillating function. The O(e^) term also contains two oscillatory functions. These average out to zero, and the first 
non trivial contribution to the real part of the effective action comes from the non-oscillatory term with magnitude 
weVl6 in Eq. (83). 

Let us next consider the case in which the entire evolution up to some time t is adiabatic, the oscillator being 
started off in the instantaneous vacuum state at some time to < * in the adiabatic regime. 

The initial conditions for an instantaneous vacuum correspond to setting a = 1, 6 = and /// = ico{to). For these 
initial conditions, because of Eq. (76), even at the lowest order of approximation in the adiabatic regime B cannot be 
assumed to remain constant. To a first approximation, we now have: 



A{t) « 1 , B{t) « ^e'*^(*) (84) 
2,10 



so that 



B{t) = -\ (e(t)e2^^« - e(to)) + ^ / dt'e{t')e^''^''\ (85) 

J to 

The integral over e will be again ignored under the assumption that e{t) ^ 0. So, at the lowest order of approximation, 

-W = ||e--^« = -i(e(t)-e(to)e-^^W), 

\z{t)f = ^ [e\t) + e\to) - 2e{t)e{to) cos {2e{t))] . (86) 
This can also be derived from the equation for z: 

z + 2iujz= ^{l- z^). (87) 
2u) 

Starting with z{tQ) = and assuming that to the lowest order of approximation {1 — z^) « 1 the above equation can 
be solved to obtain the same expression for z{t) as found above. The approximate expression for z gives 

rZT^ « 1 + ^ [e\t) + e\t,) - 2e{t)e{to) cos {20{t))] + 0{e% (88) 

so that the particle number to order is given by 

(n) « ^ [e^{t) + e\to) - 2e{t)e{to) cos (20(t))] . (89) 

If e(to) = 0, then the particle number, to the lowest order in e, is just e^(f)/16. 

In order to verify the validity of our approximation, let us reconsider as an illustration, the example of constant 
adiabaticity parameter we analyzed earlier. For this case, e{t) = —v and Q{f) — (l/i')ln(l + vt). [The e dependent 
integrals are exactly zero here.] This gives the approximate particle number from Eq. (89) to order as 



(n) 



1 — cos ( - ln(l + VT^ 



(90) 



This is clearly in agreement with the approximate expression derived in Eq. (68) from the exact analytic formula. 

It is also straightforward to compute the functions cr^ and J using the approximation for z in Eq. (86), and to 0{e) 
are given by 



0- = —7- I i H — sm 



{2Q(t))\ (91) 



and 



^=-^(eW-e(io)cos(20(i))). (92) 

The expressions for and J make it clear that their late time behavior depends on the asymptotic limit of the 
frequency function; in particular, if Lo(t) ^ at large times, we have ^ oo and ^ 0, so the Wigner function 
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will peak on the q-axis (this happens for example in the constant e case for <C 1). On the other hand, if uj(t) oo 
(with e{t) —> 0), then cr^ ^ while J7 is oscillatory with an increasing amplitude at late times, representing a Wigner 
function that ends up peaking on the p-axis. In contrast, the qp correlation, which has the form 

{qp)w = Ja^ ~ -\ {e{t) - e{to) cos{2d{t))) + 0{e^) (93) 

is oscillatory in general and remains finite, so the classicality parameter, too, is bounded at late times. 

For the adiabatically evolved vacuum state, the imaginary part of the effective Lagrangian again leads to the 
standard result and gives the effective number of particles produced. As for the real part, making use of Eq. (86) one 
has, up to order e^: 



KeLeff w — -—Ini 



-^(.(i)-e(to)e-^^^(*) 



^[e^(t)-e(i)e(io)cos(2^(t))]. (94) 



Using Eq. (89), one can rewrite the above expression in terms of the mean particle number, to order e , as 

ReLeff ~ oj{n) - ^e'(io) + ^e{t)e{to) cos(26'(t)). (95) 

The first term in this expression has a simple interpretation. The energy drained from the classical system due 
to particle production, —u}{n), matches with this term. (Recall that effective Lagrangian and efi^ective Hamiltonian 
differ by a sign in this case.) So to lowest order in the adiabaticity parameter, one can indeed incorporate the back 
reaction due to the production of particles using the first term in ReLg/j. The other two terms will be subdominant 
when e{to) is sufficiently small and will represent transient effects. In this context, to lowest non vanishing order in e, 
the real part of the effective Lagrangian incorporates information about particle creation. 



C. Non-adiabatic evolution at late times 



We would like to compare the features suggested by the above analysis with the alternative scenario of adiabaticity 
being violated in the late time limit. All choices for uj{t) which monotonically fall off faster than 1/t for large t 
would pass into such a phase, and for such functions, the adiabaticity parameter grows without bound at late times: 
limj^oo e{t) = 00. 

Since our focus is on understanding the influence of particle creation on the evolution towards classicality, the case 
of increasing particle number (n) is considered. The approximate asymptotic form of the solution which has this 
feature, is given (see Appendix D for an outline of the derivation) to the lowest order by 

^ = (1 + A)e*^ (96) 

with 

A ^ ciu , ^ n - cqoj - 2u)t (97) 

in which cq,Ci are real constants. With these analytic expressions, one can determine the asymptotic form of the 
mean particle number as well as the correlation: 

and 

'2{qp)w = Ja^ ^ - ----- 99 

Since it has been assumed that U!{t) ^0 at late times, the particle number clearly diverges. Further, {qp)yv is an 
increasing function of time in the t — > oo limit. Thus, steady particle creation (i.e. growth in the particle number (n)) 
at late times is accompanied by a monotonic growth of the magnitude of the q-p correlation, and so the classicality 
parameter S ends up saturating at unity. 

For this case, the late time approximations for the imaginary and real parts of the effective Lagrangian, making 
use of Eq. (97), are given by 



(100) 
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from which it is evident that the imaginary part of the effective action blows up as rapid particle production occurs. 
In comparison, the real part of L^ff falls to zero because of our constraint on the asymptotic form of (jj{t), and 
ReAgff{t) has a finite asymptotic limit. 

Based on the features revealed by the above two cases, one can begin to deUneate a relationship between the 
production of particles (increasing (n)) and emergence of classical correlations in terms of the c;lassicality parameter. 
We would also like to see how the corresponding Wigner function behaves in these two limits. In order to proceed 
with this, as well as to set the above ideas on a firmer footing, we shall next analyze two toy examples which will 
serve to demonstrate the key aspects. We will resort to a numerical approach for tracking the time evolution since 
it suffices for our purpose, which is to compare the two alternative late time scenarios and verify the validity of our 
approximations. These examples will also serve as prototypes for more realistic models, like those appearing in field 
theory in cosmological or electromagnetic backgrounds, to be discussed in Paper II [18]. Let us begin with an example 
that provides an illustration of case III B considered above. 



Example 1: Adiabatic evolution at late times 



We choose for the time dependent frequency the form u}i{t) = Xy/1 + XH'^ (— oo <t < oo). This frequency function 
is symmetric in time t. It is physically relevant too, as it appears in the study of a complex scalar field evolving in 
the background of a constant, classical electric field (see for eg. Refs. [6, 7]). 

For this function, the adiabaticity parameter given by e{t) = Xt/{1 + X'^t^)^^^ vanishes in the asymptotic limits 
{t — > ±c»), so in- and out- vacuum states can be defined in the adiabatic regions. If one starts out in the in- vacuum 
state (defined by (n) ^ as f — oo), the non-equivalence of the two asymptotically defined vacua implies that the 
state will appear populated with quanta measured with respect to the out- vacuum at late times. 

The evolution of the system can be conveniently followed in our Schrodinger picture formalism. The solution to the 
equation for z, Eq. (51), determines the complete time evolution of the quantum state (with the initial condition z = 
set well within the early time adiabatic phase at some moment \to\ >> X~^ with to < 0). The numerical solutions for 
the mean particle number and the classicality parameter as functions of the dimensionless variable t = Xt are plotted 
in figures 6(a) and 6(b). 

As anticipated, the particle number (n) starts from zero in the distant past, and saturates at a constant value in the 
T — > oo limit. However, in the intermediate region (where the adiabaticity parameter e is appreciably non-zero), (n) is 
characterized by large osillations superimposed on a steadily increasing mean, which can be interpreted as a continual 
interplay of creation and annihilation of quanta. The amplitude of these oscillations gets progressively diminished as 
the evolution proceeds into the late time adiabatic regime, and the mean settles at a constant value. The oscillations 
in (n) conform to the results we have obtained in earlier sections using analytic approximations. 

The classicality parameter <S stays very close to zero at early times, but in the late time adiabatic regime, ends 
up oscillating about 5 = as our approximate analysis in section IIIB did suggest. As is evident from the plot, the 
time-averaged mean of this oscillatory variation is zero, but the variance of S has a non-zero finite value that stays 
nearly constant in the late time limit. One can possibly interpret this as an emergence of correlations in comparison 
with the state at early times. 

One can also directly visualize the evolution of z in the complex plane. The plot for the complex trajectory of 
z is shown in figure 6(c). z starts from zero and after meandering around for a while in the intermediate phase, 
ultimately ends up, for large r, on the trajectory \z\ = constant, corresponding to a constant mean particle number. 
The oscillation of the imaginary part of z, as it circles around, is reflected in the late-time variation of the classicality 
parameter; this can be understood from the relation between z and S, Eq. (49). 

Let us look at the spreading of the Wigner function in the asymptotic limits. As is clear from the plots for and 
J in figure 7, the Wigner function starts uncorrelated and sharply peaked on the vertical p axis, and again at late 
times as well, ends up getting concentrated on the p axis. This may be compared with the corresponding classical 
trajectory, which, in both the asymptotic regions, has the approximate form 
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It is clear that in both limits, q ^ and the trajectory ends up on the p axis, coinciding with the Wigner function's 
behavior. The tracking of the classical trajectory by the Wigner function, thus, is not quite the same as the evolution 
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(a) (b) (c) 



FIG. 6: Plots of the particle number (n), classicality parameter 5 and the complex trajectory of z with time for aji(t) (e ^ 
at late times), (n) saturates at a finite value at late times, and S ends up oscillating about zero with a finite amplitude. The 
variation of S reflects the complex trajectory of z, which at late times is a circle centered on the origin of the complex plane. 




(a) (b) 

FIG. 7: Variation of the functions cr^ and J for LO\ify. The Wigner function starts uncorrelated and peaked on the p— axis, 
and again at late times ends up peaking on the p— axis. 

of the classicality parameter, which, in contrast, displays an asymmetry in time and more closely tracks the particle 
number. This example also demonstrates that interpreting classicality merely in terms of peaking on the classical 
trajectory may be suspect, since this can happen even when the oscillator is in a near-vacuum state [which is the 
case at early times] . But the variation of S suggests that the state gets appreciably more correlated in the far future 
[though remaining oscillatorily zero] in relation to early times [when it stays very nearly zero] . 

For the sake of verification of the numerical work, we have repeated the above analysis for another frequency, given 
by = X[a + b tanh(At)] (a > & > 0) which has similar early/late time limits. The analysis of this model proceeds 
exactly as in the above example. The time evolution is almost identical too, particularly with regard to the particle 
number and the q-p correlation, as is clear from the plots in figure 8. In this case, too, the mean particle number 
exhibits oscillations that get suppressed as the evolution turns adiabatic for large r. This shows that our conclusions 
are fairly generic. 

We now move on to consider a third toy example that is expected to illustrate the alternative scenario, that of 
strongly non-adiabatic evolution in the late time limit. 

Example 2: Non-adiabatic evolution at late times 

Consider the time dependent frequency having the form u!3{t) = A/(l + A^<^). This frequency function presents a 
sharp contrast to our previous examples. For this case, the adiabaticity parameter e(t) — —2Xt is small for \t\ < X~^, 
allowing a reasonable definition of an adiabatic vacuum in the region around t ~ 0. But the magnitude of e continues 
to grow linearly with time, and for large t, an adiabatic vacuum clearly cannot be defined. Thus, although the usual 
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FIG. 8: Evolution of the particle number (n), classicality parameter <S, complex trajectory of z and the Wigner function 
variables cP' and J with time for the frequency function 0)2 (t) = A[a + 6 tanh(At)] (e — > at late times) for the choice a — 2 
and 6=1. The variation of (n), S and z are qualitatively similar to those for the frequency function u)\(t') considered earlier. 

definition of particle based on an asymptotically adiabatic out- vacuum is not possible here, our choice of the particle 
number (n) provides a very reasonable alternative to quantify the state's particle content in such a situation. 

We set the oscillator in the instantaneous ground state at t = and track its evolution for i > 0. The plots 
corresponding to this solution for z are depicted in figures 9(a)-9(c). Our approximations given by Eqs. (98) and (99) 
are borne out, at least qualitatively, by the numerical plots, which show that the particle number (n) as well as the 
phase space correlation (gp)w grow rapidly without bound at late times. The state can, thus, be regarded as being 
driven to become strongly correlated with particle creation. This is also in contrast to the [asymptotically] adiabatic 
evolution for the frequency function LO\(i) [or of ^^^2(^)1 which was taken up earlier, where S stays oscillatorily zero in 
the absence of particle creation at late times. 




FIG. 9: Variation of (n), the classicality parameter S and z with time for ^jJzii) (e ^ 00 as i — > 00). (n) monotonically grows 
with time and is unbounded; S also grows sharply with increasing particle number and saturates at unity at late times. The 
trajectory of z is confined to a quadrant, and the state ends up at the point z = — 1 in the t ^ oq limit. 

The variation of z in the complex plane, shown in figure 9(c), brings out the contrast with the previous case even 
more vividly. Now z starts from the origin, stays within a quadrant, and ends up going to the limiting value of z = — 1. 
The imaginary part of z does not show oscillations, but rises to a maximum before falling off in the large t limit. 
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FIG. 10: Variation of the functions and J for LL}-i{t). The Wigner function starts uncorrelated in q and p but spread out, 
and ends up peaking on the axis at late times. 
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FIG. 11: Variation of the functions {n)uj and 5/r for the non-adiabatic case. These quantities saturate at constant values in 
the t oo limit, in agreement with our late time analytic approximations in Eq. (98) and Eq. (99). 



Let US also look at the spreading of the Wigner function at late times. From the plots for cr^ and J in figure 10, it 
is clear that the Wigner function, starting from being uncorrelated but fairly spread out in phase space at i = 0, ends 
up sharply peaking on the q axis as i — > oo. The classical trajectory, in comparison, is given in terms of the exact 
analytical solution to Eq. (16), which is expressible in terms of standard functions [24]: 
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Thus, the classical trajectory too ends up on the q axis, being tracked by the behavior of the Wigner function. In this 
instance, the interpretations of the state approaching classicality based on both peaking on the classical trajectory 
and the late time growth of {qp)w clearly match. 

The plots in fig. 9 follow our expectations based on the approximate analysis carried out in section IIIC at least 
at the qualitative level, but we would also like to check if our late time analytic approximations given by Eqs.(98) 
and (99) hold. From the expressions for (n) and S, it can be seen that at late times {t 00), the variables (n)w and 
Jq- jt approach constant limiting values. We have plotted the numerical solutions for these quantities as functions of 
time (t) in figures 11(a) and 11(b). Sure enough, they reach steady-state values for large r, in clear agreement with 
our asymptotic analysis. 
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IV. DISCUSSION 

In this paper, we have outUned a general formaUsm to analyze the dynamics of a time dependent oscillator, which 
provides a reasonable means of quantifying the physical content of the evolving quantum state, and addressing the 
issue of quantum-to-classical transitions. Among other things, our definition of the time-dependent particle number 
(n) based on instantaneous eigenstates is a very reasonable one to adopt in the absence of adiabatically definable in 
and out vacua. (The instantaneous particle concept has been considered in field theory in the literature before [1], 
and some 'difficulties' (the prediction of vastly more particle creation than expected on physical grounds in certain 
situations etc.) have been pointed out; nevertheless, we believe that in the present context of an oscillator with general 
time dependence, its suitability outweighs its disadvantages.) An interesting feature suggested by our approximate 
analysis, in particular of adiabatic evolution, is the possible incorporation of information about particle production in 
the real part of the effective Lagrangian, having implications for the issue of studying back-reaction in semiclassical 
theory, which normally involves considering only the real part of Leff. 

We next took up a detailed study of various illustrative toy examples, in an attempt to understand more clearly the 
connection between particle creation and the approach to classicality based on the Wigner function. We take stock by 
summarizing the main results of our analysis (in particular the late time behavior of various quantities constructed 
out of the quantum state, starting from an instantaneous vacuum at the initial moment) of the different examples in 
tabular form below: 





Constant adiabatic 
parameter |e| <C 1 


Constant adiabatic 
parameter |e| ^ 1 


Adiabatic at late times 
(e(i) — * as t — oo) 


Adiabatic throughout 
(starting from vacuum) 


Non-adiabatic at late 
times (e(t) ^ oo as t ^ oo) 








A\/1 + A2t2 


|€(t)| <^ 1 at all times 


X 


i+Lt ^) 




(l+A^t^) 


in) 


Oscillatory at all times 
with constant amplitude 


Monotonically 
increases with time 


Saturates at a finite 
value asymptotically 


Remains finite and of 
order at all times 


Monotonically increases 
with time 


z(f) 


A circular trajectory in 
the complex plane with 
Imz > at all times 


Stays within a quadrant 
and ends up at a point 
with |2;| = 1 


Ends up on a circular 
trajectory centered on 
the origin 


Remains confined to a 
finite region around the 
origin for general e(t) 


Stays within a quadrant 
and ends up at the 
point z = — 1 


S{t) 


Oscillatory at all times 
with nearly constant 
finite amplitude 


Grows and saturates 
at unity at late times 


Oscillatorily zero with 
finite amplitude at 


Remains finite and of 
order e at all times 


Grows and saturates 
at unity at late times 


W{q,p) 


Starts spread out 
and ends up peaking 
on the q-axis 


Starts spread out and 
ends up peaking on 
the g-axis 


Starts peaked on the p- 
axis and again ends up 
peaking on the axis 


Ends up peaking on cither 
q- or p-axis depending on 
late time limit of u'(t) 


Starts spread out and 
ends up peaking on the 
g-axis 



We began by examining an oscillator with constant adiabaticity parameter, a model that is capable of providing 
a prototype for adiabatic as well as non-adiabatic evolution. Going from small to large values of |e|, there is a 
clear variation in the features exhibited by the evolving quantum state. In the adiabatic regime, particle creation is 
suppressed and the mean particle number, starting from an instantaneous vacuum, maintains an oscillatory profile. 
The classicality parameter, taken as a measure of phase space correlations, also starts from zero (corresponding to 
the instantaneous vacuum) and remains bounded and oscillatory. This behavior is sharply contrasted by that in the 
extreme non-adiabatic case, where the mean number of particles as well as the q-p correlation diverge at late times 
(although the energy saturates at a finite limiting value). We subsequently analyzed three toy models that provide 
examples of frequency functions which vary adiabatically (|e| 0) or non-adiabatically (|e| ^ 1) at late times. These 
models provide adequate confirmation of the validity and generality of the features revealed by the constant e case. 
Our examples also demonstrate that peaking of the Wigner distribution on the classical trajectory is a rather general 
feature and occurs whenever the frequency goes to zero or infinity, and is quite independent of the particle production 
(this is particularly clear in the cases of constant |e| <C 1 and the frequency function W3(i)). 

The question of interpreting classical behavior based on the Wigner function involves considering two possible 
requirements: peaking on the corresponding classical phase space trajectory and emergence of correlation between q 
and p. To repeat once again what we have stated before, our toy examples provide indications that an interpretation 
based solely on the former requirement can leads to ambiguities. (Another example of such behavior is found in the 
cosmological context to be discussed in Paper II [18]: when evolution of scalar field modes in a de Sitter universe is 
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considered, the Wigner function can be shown to get peaked on the classical trajectory not just at late times, but also 
in the asymptotic past when the modes are started off in the Bunch-Davies vacuum state.) In light of this fact, we 
believe that the correlation function S comes in as a useful additional variable to quantify classicality. The Wigner 
function for a vacuum state is expressible as an uncorrelated product of the form f{q)g{p), and for this state, S is 
indeed zero. 

Once one chooses to concentrate on the variation of the q-p correlation, two possible asymptotic limits can be 
identified. When the system evolves adiabatically, the correlation function maintains an oscillatory profile. So 
suppression of particle creation keeps the correlation from varying monotonically (and in fact it can even average 
out to zero in some cases, like in the example exhibiting late time adiabaticity) . If, on the other hand, the system 
violates adiabaticity at late times (with w, ^ as t ^ oo), large amount of particle production occurs, and this is 
accompanied by sharp growth in the classicality parameter. This is in reasonable agreement with, and to some extent 
a generalization of, previous analysis in the literature [10, 16, 23] associating the notion of classicality with particle 
creation in the context of inflationary cosmology. 

We would here like to clarify our point of view regarding the interpretation of classicality in relation to the other 
approaches found in the literature that try to explain emergence of classical properties in terms of the behavior of 
Wigner functions. It goes without saying that the real world is fundamentally quantum mechanical, and the idea 
of classical behavior is essentially a matter of interpretation; as a consequence, there is no reason to expect a single 
criterion to determine classicality to be applicable in all situations. Different criteria are expected to be appropriate 
depending upon the context. For a normal oscillator, a WKB state is an energy eigenstate, being a solution of 
the time-independent Schrodinger equation. For the WKB state, it can be shown that the corresponding Wigner 
function is peaked on the curve P = f{Q), which represents the classical phase space trajectory [15], and thus shows 
a precise correlation between position and momentum. The Wigner function in this case does not involve explicit 
time dependence, and the above-mentioned feature is also naturally independent of time. The other state which has 
been analyzed in this context is a coherent state which is explicitly time-dependent and for which is actually 
peaked about a classical trajectory. For such a state, the Wigner function is expressible as an uncorrelated product of 
gaussians of the form \V{q,p,t) = A{q,t)B{p,t) and is time-dependent, but the peak follows a classical trajectory in 
phase space [10, 15]. So, although at any given moment, the Wigner function is not peaked on a classical trajectory 
(as happens in the case of the WKB stationary state mentioned above), the peak traces out a classical trajectory over 
time. 

The gaussian state we are working with is, firstly, explicitly time dependent (and so is the corresponding Wigner 
function) and the 'peak', corresponding to the maximum value of W, remains fixed at the origin (q — p — 0), since 
we have chosen a quantum state having zero mean. It is therefore more appropriate to consider the way the Wigner 
function is concentrated in some region of phase space (by, say, tracking the behavior of a Wigner function ellipse, 
i.e. a contour corresponding to a chosen value of W). The Wigner distribution is not concentrated on the classical 
trajectory at all times in the manner a stationary WKB state is, and does not track the classical path (unlike the 
coherent state). But in some particular limit, like at very late times, one can make a clear-cut correspondence between 
the concentration of the Wigner function, based on the orientation and shape of an ellipse, and the behavior of a 
general classical phase space trajectory in that limit. (Based on a careful analysis of the graphs for the evolution of the 
Wigner function ellipses and the classical trajectory, such a correspondence does not appear to hold at intermediate 
times.) It is only in this sense that the Wigner function gets 'peaked' on the classical trajectory in some limit, and 
our interpretations in the various examples studied have been based on making this kind of correspondence. 

The essential features revealed by our analysis of the toy models are expected to hold more generally, at least at 
the qualitative level. In Paper II [18], we shall consider examples from field theory in cosmological and electric field 
backgrounds. 
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APPENDIX 

A. Derivation of the generating function for tp{q,t) 

The amplitude for the oscillator to be in the nth instantaneous eigenstate at time t is non-zero only for even 

n, since V is an even function, and is given by 

Cn = r dg4>*n{,,tmq,t) = N ("^y^' r rf,if„(Vmc.,)e-(^+-)«^+^^c("+i)-(*)'^*. (106) 

J-oo V TT / V2"n! J-oo 

This can be evaluated by making use of the following generating function for Hermite polynomials [25]: 



.f;^. ,107, 



-t^+2tx 



n=0 

,21 



Substituting x = ^/rruoq, multiplying both sides by exp[— (i? + mu)/2)(f], and then integrating w.r.t. q finally gives 

OO . 71 7- 



Equating coefficients of equal powers of t on both sides gives 

/.. = w/=f=i^=^ (n = 0,1.2,.., ,109) 

while /„ vanishes for odd n. Therefore, the probability for the oscillator to be in the eigenstate (^2n at time t is simply 

p |2-p(2^)'N'" .nil 



and 



where 



using Eq. (18), and 



P=!S^. (113) 



The generating function for this probability distribution is defined as 

G(x).f:P2„."=pf:g^^x". (114) 

n=0 ra=0 ^ 

It can be evaluated using the following relation [25] : 

^(2n)!X"_ 

n=0 ^ ' 

which gives (setting x\z\'^ = X): 

G{x) = P {1 - x\z\'^)~^^^ . (116) 
And finally, since the total probability given by G{1) is unity, we have 

P^x/T^. (117) 
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B. Some properties of the probability distribution for the gaussian state 

Consider the probability distribution of the particle number for our gaussian quantum state, given by Eq. (23). 
Setting \z\'^ = exp (—6) and x = exp (—A), the generating function in Eq. (23) can be written as 



G(A) = yi^^. (118) 

(It may be remembered that this is a generating function for pairs of particles, i.e. the fcth term in the series 

expansion in powers of exp (—A) gives the probability for k pairs.) In terms of h, the mean particle number has the 
form in) = (exp (6) — 1)~^, which is same as for a Planckian distribution specified by a temperature . It's also 
straightforward to show that 

(n^) = 4 k'P,, = =3{nf + 2{n) (119) 

f-' dX^ |A=o 

and so the dispersion in the particle number is given by 

(An)2 = (n2) - (n)2 = 2(n) {{v) + 1) . (120) 

Let us also look at a thermal state in comparison, which is given by the following distribution: 

P(n)=jVe-''" (n = 0,l,2,...) (121) 

with the normalization N = (1 — e"*") fixed by requiring the total probability to be unity. This distribution is described 
by the generating function 

G(A) = ^P(n)e-^"= — ^ . (122) 

n=0 ^ ' 

from which it is clear that the generating functions of the two distributions are related in a simple manner as follows: 



G(A) = ^G(A). (123) 
For the thermal distribution, the first two moments are given by 

{n) = , (n^) = + ^ (124) 

^ ' - 1 ^ ' - 1 (efe - 1)^ 

and hence 

(An)2 = (n) ((n) + 1) (125) 

which is the standard result [23]. The distribution in Eq. (118), thus, has a dispersion that differs by only a factor of 
2 from that for the thermal state. This fact has been pointed out in the literature before (see for e.g. [26]), but its 
significance, if any, is unclear. 

C. Connecting with the Heisenberg picture 

For completeness and for reference, we give below an outline of the analysis of the time-dependent oscillator in the 

Heisenberg picture. 

The starting point is the Hamiltonian, defined as 

Hiq,p) = ^(p2 + a;2(f)g2) ^ith p = q (126) 

and q satisfies the equation of motion 

q{t) + iv^t)q{t) = 0. (127) 
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q is elevated to the status of an operator, and we write it as 

q{t) = af{t) + a}r{t) (128) 

where a, a) are the time independent creation and annihilation operators, and f{t) satisfies the same equation as q{t). 
Further, requiring [q.,p] = i and the ladder operators to satisfy the commutation relation [a, o^] = 1, one can show 
that 

rf-fr = i. (129) 

Next, the function / is written in terms of two new functions A{t) and B{t) in the following form: 

f{t) = A{t)e-'P^*^ + B(i)e'^W (p = u) (130) 

Since two functions have been introduced, an additional constraint needs to be specified to fix the form of these 
functions. We impose the requirement that f{t) be of the form 

f{t) = -iu;A{t)e-'P^*^ + iu>B{t)e'P^*^ (131) 

i.e. the derivative of / have the same form as it would if A, B were to be independent of time. This gives the following 
constraint on A, B: 

Ae-'P + Be'P = 0. (132) 
Further, substituting the expression in Eq. (130) in the equation for / and using (131) gives 

Ae-'P - Be'P = - {Be'f - Ae'^P) (133) 

CO 

Using Eq. (132) and Eq. (133), one ends up with the following pair of coupled equations for A and B: 

B+—B= —Ae-^^P, 

A+—A = —Be^'P. (134) 



Let A{t) = a{t) /V2u; and B{t) = /3{t) /V2to. Then Eq. (134) reduce to 

a = ^pe^'f and 

2U! 



P = ^ae-^'P. (135) 

Given the functional form of oj{t), these coupled equations specify the complete time evolution of the system starting 
from any initial condition on a and p. 

Eq. (128) is next rewritten in terms of a, (3: 

qit) = am + a^nt) = f^a+^at)e-+fMa+^at)e- 

\x/2oj x/2uj J \V2co V2oj J 



^ a{t)e-'P + a+{t)e'P _ ^^^^^ 
\f2u) 

The last line above defines for us the time-dependent creation and annihilation operators as 

a{t) = a{t)a + P*{t)a^ , a^{t) = a*{t)a^ + (3{t)a (137) 

and these imply 

a = a*{t)a{t)- (i*{t)a){t) , a'< = a{t)a'< {t) - (3{t)a{t). (138) 
{a{t),a\t) axe defined in such a way that for the constant frequency case, they remain independent of time.) 
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Let stand for the n-particle state at time t. The vacmini at time t is defined by the condition a{t)\0,t) = 0. 
Let a(0) = 1 and /3(0) = so that a and coincide with the ladder operators at t = 0. One can then expand the 
vacuum state at t = (the in vacuum) in terms of the basis of states defined at t: 

\0,0) = ^Cn\n,t). (139) 

n=0 

Operating with a on both sides and using the relation in Eq. (138) gives 

Ci = ; C2m+i = (m = 0,l,2,...) (140) 

and 



^^'"=^V^^^'™-^ (m=l,2,...). (141) 



The recurrence relation above allows one to express in terms of Cq, and thus 



0'0) = ^oE(£)'"^|2-,^>- (142) 



TO=0 



Using the relation (a^(t)) \0,t) = i/(2m)!|2m,f), the expression above can be rewritten in the following form: 



m=0 

The magnitude of Co can be set by normalization: 

|^2m| — J- =^ — 

m=0 



(0,0|0,0) = Yl l^^^r = 1 ^ -r^^= = 1 ^ l^oHal = 1 (144) 







2 






a 





so |Co| = \q\-'/^. 

One can compare the coefficient in Eq. (143) with the expression for C2m derived earlier in the Schrodinger picture, 
Eq. (110), to make the correspondence ^exp(2ip) = l3*/a* with p = ui. This relation can also be deduced by working 
out the equation satisfied by {f3*/a*) exp(— 2ip) using Eq. (135), which turns out to be the same as the z equation in 
(51). 



D. Late time analytic approximation for z when e{t) — > oo 

Making the substitution z = {1 + A)exp(i^) (with A real) in the equation for z, eq.(51), the following coupled 
equations for A & are obtained: 



uje' + 



1 + 



2(A+1) 



sin 6*+ - = 0, 
e 



wA'- 



A- 



A^ 



cos 61 = 



(145) 



where the prime (') denotes differentiation w.r.t. u). We axe seeking an approximate late time solution with A ^ 
as w — > (i.e. t — > oo). Assuming that in this limit, the O(A^) terms in eq.(145) can be dropped, we are left with 



a;6»' + sin6'+ - « 0, 
e 

u)A' + A cos 6* PS 0. 



(146) 



Let sin0 ^ as w ^ 0, i.e. ^ nir {n e I). This would give cos^ ±1. Consider the upper limit first. Then 
9 mr with even n. This gives 



wA' + A « 



(147) 
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which, in the lu ^ limit, blows up, and wc run into an inconsistency with the dropping of the O(A^) terms earlier. 
Consider therefore the other limit, i.e. cos9 —1. Let = nn — 4> {n is odd). Eqs.(146) reduce to 

/ 2 

-uxj) + pa -- , 
e 

wA' - A w 0. (148) 

The second equation above can be integrated to give A ciiv which yields the right late time limit. The first equation 
can also be solved, and gives 

- « co+ / — T,doJ = co + 2t (149) 
LJ J e ijj-' 

using the relation e = w/w^. 
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